3.3.66 \(\int \sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} \, dx\) [266]

Optimal. Leaf size=123 \[ \sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x+\frac {(b c+a d) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {a+\frac {b}{x}}}{\sqrt {a} \sqrt {c+\frac {d}{x}}}\right )}{\sqrt {a} \sqrt {c}}-2 \sqrt {b} \sqrt {d} \tanh ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b} \sqrt {c+\frac {d}{x}}}\right ) \]

[Out]

(a*d+b*c)*arctanh(c^(1/2)*(a+b/x)^(1/2)/a^(1/2)/(c+d/x)^(1/2))/a^(1/2)/c^(1/2)-2*arctanh(d^(1/2)*(a+b/x)^(1/2)
/b^(1/2)/(c+d/x)^(1/2))*b^(1/2)*d^(1/2)+x*(a+b/x)^(1/2)*(c+d/x)^(1/2)

________________________________________________________________________________________

Rubi [A]
time = 0.11, antiderivative size = 123, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 8, integrand size = 23, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.348, Rules used = {382, 99, 163, 65, 223, 212, 95, 214} \begin {gather*} x \sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}}+\frac {(a d+b c) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {a+\frac {b}{x}}}{\sqrt {a} \sqrt {c+\frac {d}{x}}}\right )}{\sqrt {a} \sqrt {c}}-2 \sqrt {b} \sqrt {d} \tanh ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b} \sqrt {c+\frac {d}{x}}}\right ) \end {gather*}

Antiderivative was successfully verified.

[In]

Int[Sqrt[a + b/x]*Sqrt[c + d/x],x]

[Out]

Sqrt[a + b/x]*Sqrt[c + d/x]*x + ((b*c + a*d)*ArcTanh[(Sqrt[c]*Sqrt[a + b/x])/(Sqrt[a]*Sqrt[c + d/x])])/(Sqrt[a
]*Sqrt[c]) - 2*Sqrt[b]*Sqrt[d]*ArcTanh[(Sqrt[d]*Sqrt[a + b/x])/(Sqrt[b]*Sqrt[c + d/x])]

Rule 65

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 95

Int[(((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_))/((e_.) + (f_.)*(x_)), x_Symbol] :> With[{q = Denomin
ator[m]}, Dist[q, Subst[Int[x^(q*(m + 1) - 1)/(b*e - a*f - (d*e - c*f)*x^q), x], x, (a + b*x)^(1/q)/(c + d*x)^
(1/q)], x]] /; FreeQ[{a, b, c, d, e, f}, x] && EqQ[m + n + 1, 0] && RationalQ[n] && LtQ[-1, m, 0] && SimplerQ[
a + b*x, c + d*x]

Rule 99

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.), x_Symbol] :> Simp[(a + b*
x)^(m + 1)*(c + d*x)^n*((e + f*x)^p/(b*(m + 1))), x] - Dist[1/(b*(m + 1)), Int[(a + b*x)^(m + 1)*(c + d*x)^(n
- 1)*(e + f*x)^(p - 1)*Simp[d*e*n + c*f*p + d*f*(n + p)*x, x], x], x] /; FreeQ[{a, b, c, d, e, f}, x] && LtQ[m
, -1] && GtQ[n, 0] && GtQ[p, 0] && (IntegersQ[2*m, 2*n, 2*p] || IntegersQ[m, n + p] || IntegersQ[p, m + n])

Rule 163

Int[(((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_)*((g_.) + (h_.)*(x_)))/((a_.) + (b_.)*(x_)), x_Symbol]
:> Dist[h/b, Int[(c + d*x)^n*(e + f*x)^p, x], x] + Dist[(b*g - a*h)/b, Int[(c + d*x)^n*((e + f*x)^p/(a + b*x)
), x], x] /; FreeQ[{a, b, c, d, e, f, g, h, n, p}, x]

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 223

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Subst[Int[1/(1 - b*x^2), x], x, x/Sqrt[a + b*x^2]] /; FreeQ[{a,
b}, x] &&  !GtQ[a, 0]

Rule 382

Int[((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> -Subst[Int[(a + b/x^n)^p*((c +
d/x^n)^q/x^2), x], x, 1/x] /; FreeQ[{a, b, c, d, p, q}, x] && NeQ[b*c - a*d, 0] && ILtQ[n, 0]

Rubi steps

\begin {align*} \int \sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} \, dx &=-\text {Subst}\left (\int \frac {\sqrt {a+b x} \sqrt {c+d x}}{x^2} \, dx,x,\frac {1}{x}\right )\\ &=\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x-\text {Subst}\left (\int \frac {\frac {1}{2} (b c+a d)+b d x}{x \sqrt {a+b x} \sqrt {c+d x}} \, dx,x,\frac {1}{x}\right )\\ &=\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x-(b d) \text {Subst}\left (\int \frac {1}{\sqrt {a+b x} \sqrt {c+d x}} \, dx,x,\frac {1}{x}\right )-\frac {1}{2} (b c+a d) \text {Subst}\left (\int \frac {1}{x \sqrt {a+b x} \sqrt {c+d x}} \, dx,x,\frac {1}{x}\right )\\ &=\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x-(2 d) \text {Subst}\left (\int \frac {1}{\sqrt {c-\frac {a d}{b}+\frac {d x^2}{b}}} \, dx,x,\sqrt {a+\frac {b}{x}}\right )-(b c+a d) \text {Subst}\left (\int \frac {1}{-a+c x^2} \, dx,x,\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {c+\frac {d}{x}}}\right )\\ &=\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x+\frac {(b c+a d) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {a+\frac {b}{x}}}{\sqrt {a} \sqrt {c+\frac {d}{x}}}\right )}{\sqrt {a} \sqrt {c}}-(2 d) \text {Subst}\left (\int \frac {1}{1-\frac {d x^2}{b}} \, dx,x,\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {c+\frac {d}{x}}}\right )\\ &=\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x+\frac {(b c+a d) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {a+\frac {b}{x}}}{\sqrt {a} \sqrt {c+\frac {d}{x}}}\right )}{\sqrt {a} \sqrt {c}}-2 \sqrt {b} \sqrt {d} \tanh ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b} \sqrt {c+\frac {d}{x}}}\right )\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 0.35, size = 176, normalized size = 1.43 \begin {gather*} \frac {\sqrt {a+\frac {b}{x}} \sqrt {c+\frac {d}{x}} x \left ((b c+a d) \sqrt {b+a x} \tanh ^{-1}\left (\frac {\sqrt {a} \sqrt {d+c x}}{\sqrt {c} \sqrt {b+a x}}\right )+\sqrt {a} \sqrt {c} \left ((b+a x) \sqrt {d+c x}-2 \sqrt {b} \sqrt {d} \sqrt {b+a x} \tanh ^{-1}\left (\frac {\sqrt {b} \sqrt {d+c x}}{\sqrt {d} \sqrt {b+a x}}\right )\right )\right )}{\sqrt {a} \sqrt {c} (b+a x) \sqrt {d+c x}} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[Sqrt[a + b/x]*Sqrt[c + d/x],x]

[Out]

(Sqrt[a + b/x]*Sqrt[c + d/x]*x*((b*c + a*d)*Sqrt[b + a*x]*ArcTanh[(Sqrt[a]*Sqrt[d + c*x])/(Sqrt[c]*Sqrt[b + a*
x])] + Sqrt[a]*Sqrt[c]*((b + a*x)*Sqrt[d + c*x] - 2*Sqrt[b]*Sqrt[d]*Sqrt[b + a*x]*ArcTanh[(Sqrt[b]*Sqrt[d + c*
x])/(Sqrt[d]*Sqrt[b + a*x])])))/(Sqrt[a]*Sqrt[c]*(b + a*x)*Sqrt[d + c*x])

________________________________________________________________________________________

Maple [B] Leaf count of result is larger than twice the leaf count of optimal. \(217\) vs. \(2(95)=190\).
time = 0.05, size = 218, normalized size = 1.77

method	result	size
default	\(\frac {\sqrt {\frac {a x +b}{x}}\, x \sqrt {\frac {c x +d}{x}}\, \left (\ln \left (\frac {2 a c x +2 \sqrt {\left (c x +d \right ) \left (a x +b \right )}\, \sqrt {a c}+a d +b c}{2 \sqrt {a c}}\right ) \sqrt {b d}\, a d +\ln \left (\frac {2 a c x +2 \sqrt {\left (c x +d \right ) \left (a x +b \right )}\, \sqrt {a c}+a d +b c}{2 \sqrt {a c}}\right ) \sqrt {b d}\, b c -2 b d \ln \left (\frac {a d x +b c x +2 \sqrt {b d}\, \sqrt {\left (c x +d \right ) \left (a x +b \right )}+2 b d}{x}\right ) \sqrt {a c}+2 \sqrt {\left (c x +d \right ) \left (a x +b \right )}\, \sqrt {a c}\, \sqrt {b d}\right )}{2 \sqrt {\left (c x +d \right ) \left (a x +b \right )}\, \sqrt {a c}\, \sqrt {b d}}\)	\(218\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((c+d/x)^(1/2)*(a+1/x*b)^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/2*((a*x+b)/x)^(1/2)*x*((c*x+d)/x)^(1/2)*(ln(1/2*(2*a*c*x+2*((c*x+d)*(a*x+b))^(1/2)*(a*c)^(1/2)+a*d+b*c)/(a*c
)^(1/2))*(b*d)^(1/2)*a*d+ln(1/2*(2*a*c*x+2*((c*x+d)*(a*x+b))^(1/2)*(a*c)^(1/2)+a*d+b*c)/(a*c)^(1/2))*(b*d)^(1/
2)*b*c-2*b*d*ln((a*d*x+b*c*x+2*(b*d)^(1/2)*((c*x+d)*(a*x+b))^(1/2)+2*b*d)/x)*(a*c)^(1/2)+2*((c*x+d)*(a*x+b))^(
1/2)*(a*c)^(1/2)*(b*d)^(1/2))/((c*x+d)*(a*x+b))^(1/2)/(a*c)^(1/2)/(b*d)^(1/2)

________________________________________________________________________________________

Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d/x)^(1/2)*(a+b/x)^(1/2),x, algorithm="maxima")

[Out]

integrate(sqrt(a + b/x)*sqrt(c + d/x), x)

________________________________________________________________________________________

Fricas [A]
time = 4.41, size = 890, normalized size = 7.24 \begin {gather*} \left [\frac {4 \, a c x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + 2 \, \sqrt {b d} a c \log \left (-\frac {8 \, b^{2} d^{2} + {\left (b^{2} c^{2} + 6 \, a b c d + a^{2} d^{2}\right )} x^{2} - 4 \, {\left (2 \, b d x + {\left (b c + a d\right )} x^{2}\right )} \sqrt {b d} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + 8 \, {\left (b^{2} c d + a b d^{2}\right )} x}{x^{2}}\right ) + \sqrt {a c} {\left (b c + a d\right )} \log \left (-8 \, a^{2} c^{2} x^{2} - b^{2} c^{2} - 6 \, a b c d - a^{2} d^{2} - 4 \, {\left (2 \, a c x^{2} + {\left (b c + a d\right )} x\right )} \sqrt {a c} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} - 8 \, {\left (a b c^{2} + a^{2} c d\right )} x\right )}{4 \, a c}, \frac {4 \, a c x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + 4 \, \sqrt {-b d} a c \arctan \left (\frac {{\left (2 \, b d x + {\left (b c + a d\right )} x^{2}\right )} \sqrt {-b d} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}}}{2 \, {\left (a b c d x^{2} + b^{2} d^{2} + {\left (b^{2} c d + a b d^{2}\right )} x\right )}}\right ) + \sqrt {a c} {\left (b c + a d\right )} \log \left (-8 \, a^{2} c^{2} x^{2} - b^{2} c^{2} - 6 \, a b c d - a^{2} d^{2} - 4 \, {\left (2 \, a c x^{2} + {\left (b c + a d\right )} x\right )} \sqrt {a c} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} - 8 \, {\left (a b c^{2} + a^{2} c d\right )} x\right )}{4 \, a c}, \frac {2 \, a c x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + \sqrt {b d} a c \log \left (-\frac {8 \, b^{2} d^{2} + {\left (b^{2} c^{2} + 6 \, a b c d + a^{2} d^{2}\right )} x^{2} - 4 \, {\left (2 \, b d x + {\left (b c + a d\right )} x^{2}\right )} \sqrt {b d} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + 8 \, {\left (b^{2} c d + a b d^{2}\right )} x}{x^{2}}\right ) - \sqrt {-a c} {\left (b c + a d\right )} \arctan \left (\frac {2 \, \sqrt {-a c} x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}}}{2 \, a c x + b c + a d}\right )}{2 \, a c}, \frac {2 \, a c x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}} + 2 \, \sqrt {-b d} a c \arctan \left (\frac {{\left (2 \, b d x + {\left (b c + a d\right )} x^{2}\right )} \sqrt {-b d} \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}}}{2 \, {\left (a b c d x^{2} + b^{2} d^{2} + {\left (b^{2} c d + a b d^{2}\right )} x\right )}}\right ) - \sqrt {-a c} {\left (b c + a d\right )} \arctan \left (\frac {2 \, \sqrt {-a c} x \sqrt {\frac {a x + b}{x}} \sqrt {\frac {c x + d}{x}}}{2 \, a c x + b c + a d}\right )}{2 \, a c}\right ] \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d/x)^(1/2)*(a+b/x)^(1/2),x, algorithm="fricas")

[Out]

[1/4*(4*a*c*x*sqrt((a*x + b)/x)*sqrt((c*x + d)/x) + 2*sqrt(b*d)*a*c*log(-(8*b^2*d^2 + (b^2*c^2 + 6*a*b*c*d + a
^2*d^2)*x^2 - 4*(2*b*d*x + (b*c + a*d)*x^2)*sqrt(b*d)*sqrt((a*x + b)/x)*sqrt((c*x + d)/x) + 8*(b^2*c*d + a*b*d
^2)*x)/x^2) + sqrt(a*c)*(b*c + a*d)*log(-8*a^2*c^2*x^2 - b^2*c^2 - 6*a*b*c*d - a^2*d^2 - 4*(2*a*c*x^2 + (b*c +
a*d)*x)*sqrt(a*c)*sqrt((a*x + b)/x)*sqrt((c*x + d)/x) - 8*(a*b*c^2 + a^2*c*d)*x))/(a*c), 1/4*(4*a*c*x*sqrt((a
*x + b)/x)*sqrt((c*x + d)/x) + 4*sqrt(-b*d)*a*c*arctan(1/2*(2*b*d*x + (b*c + a*d)*x^2)*sqrt(-b*d)*sqrt((a*x +
b)/x)*sqrt((c*x + d)/x)/(a*b*c*d*x^2 + b^2*d^2 + (b^2*c*d + a*b*d^2)*x)) + sqrt(a*c)*(b*c + a*d)*log(-8*a^2*c^
2*x^2 - b^2*c^2 - 6*a*b*c*d - a^2*d^2 - 4*(2*a*c*x^2 + (b*c + a*d)*x)*sqrt(a*c)*sqrt((a*x + b)/x)*sqrt((c*x +
d)/x) - 8*(a*b*c^2 + a^2*c*d)*x))/(a*c), 1/2*(2*a*c*x*sqrt((a*x + b)/x)*sqrt((c*x + d)/x) + sqrt(b*d)*a*c*log(
-(8*b^2*d^2 + (b^2*c^2 + 6*a*b*c*d + a^2*d^2)*x^2 - 4*(2*b*d*x + (b*c + a*d)*x^2)*sqrt(b*d)*sqrt((a*x + b)/x)*
sqrt((c*x + d)/x) + 8*(b^2*c*d + a*b*d^2)*x)/x^2) - sqrt(-a*c)*(b*c + a*d)*arctan(2*sqrt(-a*c)*x*sqrt((a*x + b
)/x)*sqrt((c*x + d)/x)/(2*a*c*x + b*c + a*d)))/(a*c), 1/2*(2*a*c*x*sqrt((a*x + b)/x)*sqrt((c*x + d)/x) + 2*sqr
t(-b*d)*a*c*arctan(1/2*(2*b*d*x + (b*c + a*d)*x^2)*sqrt(-b*d)*sqrt((a*x + b)/x)*sqrt((c*x + d)/x)/(a*b*c*d*x^2
+ b^2*d^2 + (b^2*c*d + a*b*d^2)*x)) - sqrt(-a*c)*(b*c + a*d)*arctan(2*sqrt(-a*c)*x*sqrt((a*x + b)/x)*sqrt((c*
x + d)/x)/(2*a*c*x + b*c + a*d)))/(a*c)]

________________________________________________________________________________________

Sympy [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \sqrt {a + \frac {b}{x}} \sqrt {c + \frac {d}{x}}\, dx \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d/x)**(1/2)*(a+b/x)**(1/2),x)

[Out]

Integral(sqrt(a + b/x)*sqrt(c + d/x), x)

________________________________________________________________________________________

Giac [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {could not integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((c+d/x)^(1/2)*(a+b/x)^(1/2),x, algorithm="giac")

[Out]

integrate(sqrt(a + b/x)*sqrt(c + d/x), x)

________________________________________________________________________________________

Mupad [B]
time = 22.22, size = 2500, normalized size = 20.33 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a + b/x)^(1/2)*(c + d/x)^(1/2),x)

[Out]

atan(((b*d)^(1/2)*(2*(b*d)^(1/2)*(2*(b*d)^(1/2)*(2*(b*d)^(1/2)*((2*(4*a^(9/2)*b^9*c^(19/2) - 4*a^(13/2)*b^7*c^
(15/2)*d^2 - 4*a^(15/2)*b^6*c^(13/2)*d^3 + 4*a^(19/2)*b^4*c^(9/2)*d^5))/(a^7*c^7*d^9) - (((a + b/x)^(1/2) - a^
(1/2))*(32*a^4*b^9*c^10 - 120*a^5*b^8*c^9*d + 288*a^6*b^7*c^8*d^2 - 400*a^7*b^6*c^7*d^3 + 288*a^8*b^5*c^6*d^4
- 120*a^9*b^4*c^5*d^5 + 32*a^10*b^3*c^4*d^6))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2)))) - (2*(8*a^5*b^9*c^9
*d + 16*a^6*b^8*c^8*d^2 - 48*a^7*b^7*c^7*d^3 + 16*a^8*b^6*c^6*d^4 + 8*a^9*b^5*c^5*d^5))/(a^7*c^7*d^9) + (((a +
b/x)^(1/2) - a^(1/2))*(16*a^(7/2)*b^10*c^(21/2) - 76*a^(9/2)*b^9*c^(19/2)*d + 228*a^(11/2)*b^8*c^(17/2)*d^2 -
168*a^(13/2)*b^7*c^(15/2)*d^3 - 168*a^(15/2)*b^6*c^(13/2)*d^4 + 228*a^(17/2)*b^5*c^(11/2)*d^5 - 76*a^(19/2)*b
^4*c^(9/2)*d^6 + 16*a^(21/2)*b^3*c^(7/2)*d^7))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2)))) - (2*(a^(7/2)*b^11
*c^(21/2) + 16*a^(9/2)*b^10*c^(19/2)*d - 42*a^(11/2)*b^9*c^(17/2)*d^2 + 25*a^(13/2)*b^8*c^(15/2)*d^3 + 25*a^(1
5/2)*b^7*c^(13/2)*d^4 - 42*a^(17/2)*b^6*c^(11/2)*d^5 + 16*a^(19/2)*b^5*c^(9/2)*d^6 + a^(21/2)*b^4*c^(7/2)*d^7)
)/(a^7*c^7*d^9) + (((a + b/x)^(1/2) - a^(1/2))*(146*a^4*b^10*c^10*d - 556*a^5*b^9*c^9*d^2 + 1006*a^6*b^8*c^8*d
^3 - 1192*a^7*b^7*c^7*d^4 + 1006*a^8*b^6*c^6*d^5 - 556*a^9*b^5*c^5*d^6 + 146*a^10*b^4*c^4*d^7))/(2*a^7*c^7*d^9
*((c + d/x)^(1/2) - c^(1/2)))) + (2*(2*a^4*b^11*c^10*d + 8*a^5*b^10*c^9*d^2 - 2*a^6*b^9*c^8*d^3 - 16*a^7*b^8*c
^7*d^4 - 2*a^8*b^7*c^6*d^5 + 8*a^9*b^6*c^5*d^6 + 2*a^10*b^5*c^4*d^7))/(a^7*c^7*d^9) - (((a + b/x)^(1/2) - a^(1
/2))*(65*a^(7/2)*b^11*c^(21/2)*d - 297*a^(9/2)*b^10*c^(19/2)*d^2 + 597*a^(11/2)*b^9*c^(17/2)*d^3 - 365*a^(13/2
)*b^8*c^(15/2)*d^4 - 365*a^(15/2)*b^7*c^(13/2)*d^5 + 597*a^(17/2)*b^6*c^(11/2)*d^6 - 297*a^(19/2)*b^5*c^(9/2)*
d^7 + 65*a^(21/2)*b^4*c^(7/2)*d^8))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2))))*1i - (b*d)^(1/2)*(2*(b*d)^(1/
2)*(2*(b*d)^(1/2)*(2*(b*d)^(1/2)*((2*(4*a^(9/2)*b^9*c^(19/2) - 4*a^(13/2)*b^7*c^(15/2)*d^2 - 4*a^(15/2)*b^6*c^
(13/2)*d^3 + 4*a^(19/2)*b^4*c^(9/2)*d^5))/(a^7*c^7*d^9) - (((a + b/x)^(1/2) - a^(1/2))*(32*a^4*b^9*c^10 - 120*
a^5*b^8*c^9*d + 288*a^6*b^7*c^8*d^2 - 400*a^7*b^6*c^7*d^3 + 288*a^8*b^5*c^6*d^4 - 120*a^9*b^4*c^5*d^5 + 32*a^1
0*b^3*c^4*d^6))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2)))) + (2*(8*a^5*b^9*c^9*d + 16*a^6*b^8*c^8*d^2 - 48*a
^7*b^7*c^7*d^3 + 16*a^8*b^6*c^6*d^4 + 8*a^9*b^5*c^5*d^5))/(a^7*c^7*d^9) - (((a + b/x)^(1/2) - a^(1/2))*(16*a^(
7/2)*b^10*c^(21/2) - 76*a^(9/2)*b^9*c^(19/2)*d + 228*a^(11/2)*b^8*c^(17/2)*d^2 - 168*a^(13/2)*b^7*c^(15/2)*d^3
- 168*a^(15/2)*b^6*c^(13/2)*d^4 + 228*a^(17/2)*b^5*c^(11/2)*d^5 - 76*a^(19/2)*b^4*c^(9/2)*d^6 + 16*a^(21/2)*b
^3*c^(7/2)*d^7))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2)))) - (2*(a^(7/2)*b^11*c^(21/2) + 16*a^(9/2)*b^10*c^
(19/2)*d - 42*a^(11/2)*b^9*c^(17/2)*d^2 + 25*a^(13/2)*b^8*c^(15/2)*d^3 + 25*a^(15/2)*b^7*c^(13/2)*d^4 - 42*a^(
17/2)*b^6*c^(11/2)*d^5 + 16*a^(19/2)*b^5*c^(9/2)*d^6 + a^(21/2)*b^4*c^(7/2)*d^7))/(a^7*c^7*d^9) + (((a + b/x)^
(1/2) - a^(1/2))*(146*a^4*b^10*c^10*d - 556*a^5*b^9*c^9*d^2 + 1006*a^6*b^8*c^8*d^3 - 1192*a^7*b^7*c^7*d^4 + 10
06*a^8*b^6*c^6*d^5 - 556*a^9*b^5*c^5*d^6 + 146*a^10*b^4*c^4*d^7))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2))))
- (2*(2*a^4*b^11*c^10*d + 8*a^5*b^10*c^9*d^2 - 2*a^6*b^9*c^8*d^3 - 16*a^7*b^8*c^7*d^4 - 2*a^8*b^7*c^6*d^5 + 8
*a^9*b^6*c^5*d^6 + 2*a^10*b^5*c^4*d^7))/(a^7*c^7*d^9) + (((a + b/x)^(1/2) - a^(1/2))*(65*a^(7/2)*b^11*c^(21/2)
*d - 297*a^(9/2)*b^10*c^(19/2)*d^2 + 597*a^(11/2)*b^9*c^(17/2)*d^3 - 365*a^(13/2)*b^8*c^(15/2)*d^4 - 365*a^(15
/2)*b^7*c^(13/2)*d^5 + 597*a^(17/2)*b^6*c^(11/2)*d^6 - 297*a^(19/2)*b^5*c^(9/2)*d^7 + 65*a^(21/2)*b^4*c^(7/2)*
d^8))/(2*a^7*c^7*d^9*((c + d/x)^(1/2) - c^(1/2))))*1i)/((b*d)^(1/2)*(2*(b*d)^(1/2)*(2*(b*d)^(1/2)*(2*(b*d)^(1/
2)*((2*(4*a^(9/2)*b^9*c^(19/2) - 4*a^(13/2)*b^7*c^(15/2)*d^2 - 4*a^(15/2)*b^6*c^(13/2)*d^3 + 4*a^(19/2)*b^4*c^
(9/2)*d^5))/(a^7*c^7*d^9) - (((a + b/x)^(1/2) - a^(1/2))*(32*a^4*b^9*c^10 - 120*a^5*b^8*c^9*d + 288*a^6*b^7*c^
8*d^2 - 400*a^7*b^6*c^7*d^3 + 288*a^8*b^5*c^6*d^4 - 120*a^9*b^4*c^5*d^5 + 32*a^10*b^3*c^4*d^6))/(2*a^7*c^7*d^9
*((c + d/x)^(1/2) - c^(1/2)))) - (2*(8*a^5*b^9*c^9*d + 16*a^6*b^8*c^8*d^2 - 48*a^7*b^7*c^7*d^3 + 16*a^8*b^6*c^
6*d^4 + 8*a^9*b^5*c^5*d^5))/(a^7*c^7*d^9) + (((a + b/x)^(1/2) - a^(1/2))*(16*a^(7/2)*b^10*c^(21/2) - 76*a^(9/2
)*b^9*c^(19/2)*d + 228*a^(11/2)*b^8*c^(17/2)*d^2 - 168*a^(13/2)*b^7*c^(15/2)*d^3 - 168*a^(15/2)*b^6*c^(13/2)*d
^4 + 228*a^(17/2)*b^5*c^(11/2)*d^5 - 76*a^(19/2)*b^4*c^(9/2)*d^6 + 16*a^(21/2)*b^3*c^(7/2)*d^7))/(2*a^7*c^7*d^
9*((c + d/x)^(1/2) - c^(1/2)))) - (2*(a^(7/2)*b^11*c^(21/2) + 16*a^(9/2)*b^10*c^(19/2)*d - 42*a^(11/2)*b^9*c^(
17/2)*d^2 + 25*a^(13/2)*b^8*c^(15/2)*d^3 + 25*a^(15/2)*b^7*c^(13/2)*d^4 - 42*a^(17/2)*b^6*c^(11/2)*d^5 + 16*a^
(19/2)*b^5*c^(9/2)*d^6 + a^(21/2)*b^4*c^(7/2)*d^7))/(a^7*c^7*d^9) + (((a + b/x)^(1/2) - a^(1/2))*(146*a^4*b^10
*c^10*d - 556*a^5*b^9*c^9*d^2 + 1006*a^6*b^8*c^8*d^3 - 1192*a^7*b^7*c^7*d^4 + 1006*a^8*b^6*c^6*d^5 - 556*a^9*b
^5*c^5*d^6 + 146*a^10*b^4*c^4*d^7))/(2*a^7*c^7*...

________________________________________________________________________________________